Polarization dependence of the two-photon Pranz-Keldysh effect 
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The effect of a constant electric field on two-photon absorption in a direct band gap semiconductor 
is calculated using an independent-particle theory. Two band structure models for GaAs are used: 
a two-band parabolic model and an eight-band k - p model. Both predict a strong dependence of the 
two-photon electroabsorption spectrum on the polarization of the light with respect to the constant 
field. We attribute the polarization dependence to the strong effect of a constant field on intraband 
dynamics. 



Multiphoton absorption is an important topic in the 
study of light-matter interaction, and will only become 
more important with our increasing ability to control 
the spectral and quantum statistical properties of light. 
Pulse shaping,^ use of entangled photonspland quantum 
interference between absorption pathways^ have all been 
found to affect multiphoton absorption processes in pro- 
found and useful ways. Nonlinear absorption has also 
become an increasingly important issue in photonics!^ 
Despite strong interest in this area, the effects of constant 
electric fields have not received much attention, which is 
surprising considering that strong fields are often present 
in semiconductor-based photonic devices. Field-induced 
changes in the optical absorption in a bulk semiconduc- 
tor are predominantly caused by the Franz-Keldysh effect 
(FKE) J*"™ related to the acceleration of carriers by a con- 
stant (dc) electric field. In the case of linear absorption, 
the field induces an exponential tail for photon energies 
below the band gap caused by field-induced tunneling, 
and oscillations for energies above the band gap, related 
to the coherence of carriers. 

The two-photon FKE was the subject of a number of 
theoretical studies decades agoPKS There has recently 
been a resurgence of interest, with theoretical wor ks on 
bulk semiconductors by Garcia and Kal yanar amarP^^ 
and in nanostructures by Xia and SpectorP^^ Both used 
two-band models in the parabolic band approximation 
(PBA) to find analytical expressions for the two-photon 
absorption coefficient in the presence of an electric field. 
Experimental results for photon energies below the band 
gap have also recently been published.^ Here we extend 
a recently described theory^ of the one-photon FKE in 
bulk semiconductors to two-photon absorption. We cal- 
culate two-photon electroabsorption (EA), the difference 
between the absorption spectrum in the presence and 
in the absence of a dc field, using an eight-band k ■ p 
model for GaAs and find that the results agree qualita- 
tively with analytical expressions from a two-band PBA 
model. The theory predicts a strong dependence of the 
two-photon EA spectrum on the polarization of the opti- 
cal field with respect to the dc field. We give expressions 
for two-photon absorption with optical pulses, enabling 
calculation of effects that depend on the pulse shape, - 



and simpler expressions suitable for long pulses. 

We consider a semiconductor in the presence of a uni- 
form, constant electric field Ed c - Much of the theory is 
the same as for one-photon absorption^ here we only 
sketch the approach and give the relevant equations. In 
the independent-particle approximation, neglecting the 
interaction between carriers and scattering processes, the 
Hamiltonian including the effect of the dc field is 



n dc (t) = 



i 

2m 



-V + fiK(t) +V(x), 



where K(t) = eEd c i/S- describes the acceleration of elec- 
trons by the electric field and V(x) describes the inter- 
action of electrons with the ions of the crystal. We use 
insta ntaneo us eigenstates {</>„(k;x)} of rldcit) as basis 
statesj^2l23 denoting velocity matrix elements V m „ (k; t) , 
electric dipole matrix elements /2 m „(k;£), and band en- 
ergies w„(k). We denote destruction operators for these 
states as 6„k- To account for coupling between bands 
due to the dc field, we calculate absorption using states 
{x„(k;x)} that are related to {0 n (k;x)} by a unitary 
evolution matrix L(k; t) satisfying 



ifl 



dh (k; t) 
dt 



where T pm (k; t) 



[T(k;i)+S(k;t)]L(k;t), 
rntiLu m (k + K(t)) and S pm (k;t) 



Mpm(k; • Edc- 

It is assumed that the nominal dc field is turned on 
before the optical pulse arrives. Optical absorption is 
calculated in an interaction picture, where the ket ^(i)) 
evolves according to the interaction Hamiltonian /f e g(t), 

.<* !*(*)) 

dt 

Here H e g(t) describes the effect of the optical field 
E op t(i) = — (l/c)dA ovt /dt in the presence of the dc field; 
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it is given by 
H eS (t) 



6 Lk 6 »ik Aopt(t)- V„ 2ni (k;f) 
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where 



V nq (k; t) = J2 J4nte *)V mp (k; t)L pq (k; t). (2) 
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The initial condition of Eq. |lj) \^ H ) is the state of the 
system at the moment the dc field is turned onPSl The 
matrix elements V„ g (k;f) appearing in H c g(t) are effec- 
tive matrix elements that take into account the fact that 
the dc field can induce transitions between these bands 
itself. We assume an optical pulse of the form 



E op t (t) 



diu 
2^ 



E(w)e" 



(3) 



and assume that E(oj) does not contain frequency com- 
ponents above the band gap so one-photon absorption 
may be neglected. 

We calculate JvE^ 2 )) = ^^(oo)) in a perturbative so- 
lution of the Schrodinger equation [Eq. Q], i.e. \&(t)) = 

j^fff) + |^f(l)(t)) + |^(2)( t )) + ... ; which is 



*( 2 )\ = 



/OO pt' 
#eff(0 / H cS (t")\^ H )dt"dt'. 
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The number of carriers injected in some volume is AN = 
^\Ef( 2 )|\]>( 2 )y For the carrier density injected An, we find 

An = J duj a J duj d dJ d n l3lm (uj a , uj dl uj' d ) 



x E l uj, 



-u) d j E 3 (uj a + \^uj d 
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col, E r * 
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where i,j,l,m denote vector components and we have 
defined a carrier injection tensor 



e \ - f dk 

= 2^J 4^ 2 

cv 

in which 



^LK.^fe^,^)]*, (5) 



F c l n (kj_; -(J a + ^d)Fl v (k±; -ui a - \uj d ) 



+ 



where, finally, 

F mn (k x ;-uj) = / V mn (k x ;i)e iw *cft. 



(6) 



(7) 



Equation Q is relevant for a pulsed optical field, 
including effects due to pulse shaping.^ In the limit 
of a long optical pulse with a constant optical field 
envelope E, one can derive a Fermi golden rule 
(FGR) expression for the rate of carrier injection n = 



16w 3 r] ijlm (uj,0,0)E i E 3 [E l E m ]*. The third-order nonlin- 
ear susceptibility is given by lm[x^(— uj; — uj, uj, uj)} — 
(167r 3 ?i/3)77(w, 0, 0). For the dc field pointing along a 
crystal direction, the nonlinear absorption coefficient 
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= 21 \ 2 »7 (m.0,0), 
n z (uj)c z 
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(no summation implied) for light linearly polarized in the 
i direction. 

For two parabolic bands and a k-independent inter- 
band dipole matrix element V cl) , the equations can be 
solved analytically, and the two-photon absorption spec- 
trum may be expressed in terms of Airy functions. We 
assume that the bands have the form huj c (k) = huj g + 
Ti 2 k 2 / (2m c ) and Huj v (k) = —h 2 k 2 /(2mh), where m c is 
the electron mass, nih is the hole mass, and Tiujg is the 
band gap. The energy difference between the valence and 
conduction bands is huj cv (k) = huj g + h 2 k 2 /(2^i), where 
[i = (m^ 1 + ) _1 is the reduced mass. We assume a 
dc field pointing in the z direction and explicitly denote 
perpendicular and parallel components of the wavevec- 
tor, k = kj_ + fciiz. The off-diagonal components are 



F CT (k_L; -uj) = V OT ^Ai 



uj g — hkj_/2fj, 

n 



(9) 



where Ai(:r) is the Airy function and the electro-optic 
frequency O = (he 2 /2/i) 1 ^ 3 , defining a normalized dc field 
e ee (e/h)E dc . Using V nn (k) = hk/m n and V nn (k;t) = 
V nn (k±;t + fc||/e)pl we have *V nn (k±;t) = hk±/m n + 
zhet/m n . Using this in Eq. |2]) and Eq. ^ and using the 
fact that S = in the absence of interband coupling,^! 
we find 



F nrl (kj 



hk 



±5(uj) 



z — S (uj), 



(10) 



where we have used te~~ tult — id(e~ wt ) / du . We see that 
the diagonal (intraband) matrix element is qualitatively 
different in the direction of the dc electric field z, not 
surprising considering that the electric field accelerates 
carriers. 

We find, using Eqs. ^ and ([To]) in Eq. 
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and 



(uj a ,uj d ) = 
—8ie 2 k x uj n 



irhQ[i(4:UJ a 
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The y-component is simply the above with x — > y. 

Using Eqs. ^ and @ in Eq. Q we find in the FGR 
limit, for the optical field parallel to the dc field, 



: (^o,o) = \v*\ 



1 
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uj n — 2uj\ f uj 
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whereas for the optical field perpendicular to the dc field 
we have 



ri xxxx (u),Q,0) = 



W' 1 




-Ai 



UJg — 2CJ^ , I UJg — 2UJ 



n 



(14) 



In the limit of no dc field, one can derive using asymptotic 
expressions for Ai(x) that Eqs. (Jl3[ ) and ( [14] ) both reduce 
to the well-known expression^ for allowed-forbidden 
transitions between parabolic bands 



rT>,0,0) 



2A^ti > ' 2 u 
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-,„ A 3 / 2 
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The two-photon absorption coefficient for the two po- 
larization configurations is plotted in Fig. [[} Results 
from the PBA model are shown in Fig. [l£i, assuming 
hV* v = 10.3 eV-A and the light hole reduced mass in 
GaAs. We have multiplied the PBA results by 2 to ac- 
count for the spin degeneracy not included in the model. 

Note that the expression for fields perpendicular 
(ij 1 ™ 1 ) is the same as the one derived by Garcia,^ but 
the one for fields parallel (rj zzzz ) is very different. This 
difference is especially obvious in the two-photon elec- 
troabsorption spectrum, shown in Fig. [T}d. At a fun- 
damental level this is not surprising: Although the 
FKE goes beyond perturbation theory, consider for the 
moment a simple picture where the rfi arise from a 
perturbation due to the dc field of the imaginary part 
of X(3) (— w; —uj, uj, uj) that describes two-photon absorp- 
tion. For light polarized along the x and z directions 
the appropriate components are x[z) X an d X(|f z re ~ 
spectively, which are obviously equal for a cubic crys- 
tal such as GaAs. Since the band structure models we 
use here neglect the lack of inversion symmetry in the 
crystal, the leading order correction of these x%) due 
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FIG. 1: (color online) Calculated two-photon Franz-Keldysh 
effect in GaAs for Edc = 67 kV/cm. Blue dotted: no dc 
field; red dashed: E opt _L Ed c ; black solid: E opt || Ed c - (a) 
Two-photon absorption spectrum /3(oj) using a model with 
two bands in the parabolic band approximation (PBA, thin 
lines) and an eight-band k ■ p model (thick lines). The inset 
shows a close-up view of the absorption spectrum near the half 
bandgap Huj 9 /2, with the units on the two axes the same as 
the main plot. (b,c) Two-photon electroabsorption spectrum 
[/3(uj) with a dc field minus /3(u) without a field]. The half 
band gap is marked with an arrow, (b) Two-band PBA model, 
(c) Eight-band k ■ p model. 



to a dc field involve two powers of the field. If we let 
Yijkimn re f er f ornia lly to taking the second derivative of 

X^ 1 (—uj\—uj,uj,uj) with respect to an applied dc field, 
with components in the m and n direction, then since we 
take our electric field in the z direction, in a perturba- 



tive description we could identify T 



and 



^zzzzzz 



■q zzzz . Even within this simple picture, be- 
cause of the high rank of T we would have no reason to 
expect Y xxxxzz and r zzzzzz (and thus rj xxxx and r\ zzzz ) 
to be equal, especially since the frequencies associated 
with the last two components of the T's are drastically 
different than those of the others. 

Two-photon absorption in GaAs for energies near the 
half band gap is mostly due to two-band processes, and 
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the two-band model can be expected to capture the 
physics of that caseP^ However, three-band processes 
that involve an intermediate band also contribute to 
two-photon absorption.^ To examine how the polariza- 
tion dependence changes when intermediate bands are 
included, we calculated the two-photon EA spectrum us- 
ing a k • p modeP^ for the band structure. The model 
directly includes eight bands: the six uppermost valence 
bands and the two lowermost conduction bands. The spin 
splitting due to the lack of a center of inversion symmetry 
in GaAs is not accounted for in the eight-band model, so 
there is spin degeneracy at all k. The coupling Pq be- 
tween the valence bands and conduction bands produces 
some of the band curvature. Effects of remote bands on 
the valence band curvature are handled through modi- 
fied Luttinger parameters. A remote band parameter F 
fixes the conduction band effective mass to the observed 
value. The model parameters are the same as was used 
recently in a calculation of the one- photon FKEP The 
matrix elements and the evolution matrix are calculated 
as in the one-photon calculation. 2 ^ We show the results 
of the k • p calculation in Fig. [IJi and [TJ;. Overall, the 
results are in qualitative agreement with the two-band 
PBA model. The difference in magnitude can be traced 
to the more accurate treatment of the effective mass of 
holes in the k • p model. 

In summary, we have calculated the two-photon FKE 



and find a strong polarization dependence that we at- 
tribute to the strong effect of a dc field on intraband 
dynamics. A more sophisticated model of the band struc- 
ture, such as a 14-band modeP^ recently used for the one- 
photon FKEpSI should in addition show effects that are 
odd in the dc field because of the lack of inversion sym- 
metry in GaAs. Other future extensions to the the ory 
could include adding the Coulomb interactio n 1 11 * 13 * 1 ^ and 
scattering. Higher-order multi-photon absorption could 
be calculated by going to higher order in H e &W^ 

Previous experimental workP^ on the two-photon FKE 
only employed a dc field perpendicular to the optical 
field, and only measured the effect for Tvjj below the half 
band gap. A transverse geometry is required to observe 
the polarization dependence we predict, and that will re- 
quire application of an external field, 3 ^ rather than the 
usual use of the built-in field in a doped heterostructure. 
The prediction of a strong polarization dependence of the 
two-photon FKE, with very strong Franz-Keldysh oscilla- 
tions for all fields parallel, should stimulate more interest 
both in experiments and in more sophisticated theories 
of nonlinear absorption in external fields. 
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